Particle Condensation in Pair Exclusion Process 
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Condensation is characterized with a single macroscopic condensate whose mass is proportional 
to a system size A*'. We demonstrate how important particle interactions are in condensation phe- 
nomena. We study a modified version of the zero-range process by including a pair exclusion. Each 
particle is associated with its own partner, and particles of a pair are forbidden to stay at the 
same site. The pair exclusion is weak in that a particle interacts with only a single one among all 
others. It turns out that such a weak interaction changes the nature of condensation drastically. 
There appear a number of mesoscopic condensates: the mass of a condensate scales as mcon ~ A''^''^ 
and the number of condensates scales as A'con ~ N^^^ with a logarithmic correction. These results 
are derived analytically through a mapping to a solvable model under a certain assumption, and 
confirmed numerically. 
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Condensation occurs in various systems ranging from 
equilibrium ideal Bose gases to nonequilibrium lattice 
gases (see Ref. jl[ and references therein). In a condensed 
phase, a finite fraction of particles constitute a macro- 
scopic condensate at a single microscopic state which is 
either favored or selected spontaneously. Since particles 
are highly populated in a condensate, particle interaction 
can play an important role 0, In this paper, we in- 
vestigate robustness of a macroscopic condensate against 
a weak particle repulsion. 

Our work is motivated from a recent work on evolv- 
ing networks which exhibit a condensation of edges at 
some vertices called hubs The edge condensation 
becomes analogous to a particle condensation if one in- 
terprets vertices and edges as lattice sites and particles, 
respectively. Interestingly, the edge condensation shows 
a distinct feature that there are a large number of meso- 
scopic condensates, the origin of which has not been un- 
derstood yet. We notice that networks are subject to 
the constraint which forbids self-loops. Although it has 
a negligible effect on network properties in most cases, 
it may be responsible for the unusual condensation. We 
will investigate a lattice gas which is equipped with a so- 
called pair exclusion mimicking the self-loop constraint. 

We study the zero-range process (ZRP) Q which is 
useful for the study of condensation P, Q|. Suppose 
that there are M = pN particles in an arbitrary lat- 
tice of N sites with particle density p. A lattice site 
i = I, . . . , N can accommodate many particles, and its 
occupation number rrij will be called a mass. Each site 
i emits a particle at a mass-dependent rate u(mi). The 
canonical form is given by 



u{m) = 1 H 

m 



(1) 



with a control parameter b. An emitted particle then 
tries to hop to a site j selected with the probability given 
by a stochastic hopping matrix T = [Tji] which reflects 



the structure of the underlying lattice and the presence 
of any bias in dynamics. We now introduce the pair 
exclusion (PE): Each particle is associated with a unique 
partner so that there are M/2 distinct pairs in total with 
even M . A hopping will be rejected if it would result in 
occupation of a pair at the same site. The resulting model 
will be referred to as a pair exclusion process (PEP) in 
short. A particle interacts with its partner only. Such 
a weak interaction changes the nature of condensation 
drastically as we shall see below. 

In the ZRP without PE, the particle hopping rate de- 
pends only on the mass at departing sites. This property 
yields that the mass distribution in the stationary state is 
factorized into a product form [l| . Particle hopping rates 
in the PEP, however, depend on the whole distribution 
of particle species, which makes it difficult to treat the 
problem analytically. Hence, we adopt a series of approx- 
imations. 

We first assume that the particle species distribution is 
purely random. Then, a configuration of the system can 
be specified with a mass distribution m = (toi, . . . , m^v), 
and the particle hopping rate from a site i to j is given 

by 



Wji{mj,rn.i) = v{mj,mi) Tj^ u{mi), 



with 



{j(m, m!) = 1 — A 
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The factor v{m, to') accounts for the probability that a 
particle selected from a site of mass to' does not find its 
partner at another site of mass to. For convenience, we 
insert a parameter A = 1. One can turn off the PE by 
choosing A = 0. With < A < 1, one can deal with a pair 
repulsion that forbids a pair occupation with probability 
A. 
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For large M ^ m and m', one can expand v{m, m!) as 
{i(m, m') = v{m) + 0{{m/Mf) with 



m m 
-M = l-A^ = l-A— . 



(4) 



We make an additional approximation by ignoring the 
0{{m/MY) term. Then, the hopping rate is given by 



Wji{mj,mi) = v{mj) T.ji u{m^) 



(5) 



The hopping rate depends on the mass at both depart- 
ing and arriving sites. It is interesting to show whether 
the stationary state of such a system is also given by a 
factorized product state as the ZRP. The factorizability 
condition for the general hopping rate has been inves- 
tigated recently in an one-dimensional lattice or a com- 
plete graph [8] . The hopping rate in Eq. ([5]) is specifically 
given by the product of u(m) and v{m). We will show 
that such a system has a factorized stationary state for 
any functions u(m) and v{m) when the hopping matrix 
T satisfies the detailed balance condition. 

Let P{m; t) be the probability of the system being in 
a configuration m = (toi, . . . , m^v) at time t. The time 
evolution follows the master equation 



dP{m;t) 



N 



- {-Wjrimj,mi)P{m;t) 

+ Wj.imj -l,m, + l)P{m'j-t)} , (6) 

where m^j = {m[, . . . jTn'j^) is a shorthand notation for 
the the mass distribution with m'^ = -t- 1 , m^- = — 1 , 
and = mk for k ^ 

We seek for the stationary state solution in a product 
form given by 



Pstim) = 



1 



Z{N, M) 



N 



(7) 



where 5a,b is the Kronecker-(5 symbol, Z{N, M) is a nor- 
malization, and fi{mys have to be determined. Inserting 
Pst into the right-hand side of Eq. ([5]) and taking it to 
zero, we find that fi{m) is given by 



n=l 



;(n-l) 



u(n) 



(8) 



for m > 1 and /i(0) = 1. Here the auxiliary constant 
{uJi} should satisfy 



N 

E 
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Tij ) v{mi)u{mj) = 



For arbitrary functions u{m) and v(m), it requires 



Tji UJi — Tij ujj 



(9) 



(10) 



for all i and j. 

The stochastic matrix T by itself defines a hopping dy- 
namics of a single particle. The condition ([TU)) is equiva- 
lent to the detailed balance condition for the correspond- 
ing single particle dynamics. Hence, the PEP in one- 
dimensional rings, for example, has a factorized prod- 
uct state with {uJi = 1} when the hopping is symmetric. 
When there is a bias in hopping, the factorization is not 
guaranteed. The role of the hopping bias is discussed 
in a target process where the hopping rate depends only 
on the mass at arriving sites [8|. We note that the PEP 
in the solvable case has the same stationary state as the 
ZRP with particle emitting rate {t(m) = u{m) / v{m — 1). 

The detailed balance is a sufficient condition for the 
factorization. One may also have a factorized state with- 
out the detailed balance for a specific u{m) and v{m). 
The ZRP corresponds to the specific case with v{m) — 1. 
Using J2i'Pij = Ij can show that the relation (|9]) 
is fulfilled by {uji} satisfying J2j=i'^ij^j — ^r all 
i Q. The solution {uji} is nothing but the stationary 
state probability of a single particle whose dynamics is 
governed by T, which exists for any stochastic matrix T. 

In order to investigate the effect of the PE on the con- 
densation, we consider the PEP with the symmetric hop- 
ping in an one-dimensional ring or the random hopping 
in a complete graph. In both cases, uJt = 1 and Pst{'m) 
is given by Eq. ([7]) with fi{m) = f{m) for all i, where 



= n 



1 - A(n - l)/ipN) 



n=l 

For m ^ A^, it is given by 



1 + b/n 



r(m + i)r(5 + i) 

f(m) - — FT — , , , IX exp 
1 (to + + 1) 



2pN 



(11) 



(12) 



with the gamma function T(s) = x'^^^e^^dx. The 
single-site mass distribution function p{m), defined as the 
probability that the mass of a site is m, is given by 

p{m) = f{m)Z{N - 1,M ~m)/Z{N,M). (13) 

The PE imposes an upper cutoff at nipE — 0{N^/^) in 
the mass distribution. The knowledge on the ZRP (A = 
0) allows us to guess the phase behavior of the PEP. 
When A = 0, the system is in a fluid phase for p < 
Pc = 1/(6 — 2) and in a condensed phase for p > pc- 
The former is characterized with p{m) ~ TO~''e~™/'"" 
with a constant too = 0(1), while the latter with p{m) 
displaying a power-law scaling p(to) ^ m^^ for bulk 
sites and a peak around m — {p — Pc)N for a macro- 
scopic condensate. When the PE is turned on, the fluid 
phase will not be affected because the mass distribution 
is limited to the region to < toq where the PE is negligi- 
ble (too <C tupe)- On the other hand, the system cannot 
maintain the macroscopic condensate in the presence of 
PE because of the cutoff mpE- The excess mass cannot 
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be absorbed into bulk sites because the mean occupation 
number at bulk sites is limited hy p^. A possible solution 
is to build up (p — Pc)N/mpE mesoscopic condensates 
of mass m mpE- Hence we expect that the PE does 
not modify the phase boundary but changes the nature 
of the condensate. 

Detailed properties of the PEP can be calculated ex- 
plicitly. Interestingly, the PEP has the same stationary 
state as that of the ZRP with the hopping rate 13, U j 

h m 

uzRp{rn) 



1 



b 

m 



m 



(14) 



The ZRP has the factorized stationary state with 
fzRp{m) = ]X^^i[uzRp{n)]~^. In the large N limit, 
fzRp{fn) becomes equal to /(m) with c = A/p. The 
nature of the condensed phase is documented well in 
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Ref. [11[. Hence, we make use of the results of Ref. 
to present the property of the PEP. 

The PEP undergoes a condensation transition at the 
phase boundary p = Pc{b) = l/{b — 2). When p < pc, the 
system is in a fluid phase withy)(m) ^ m"*" exp(— m/rno)- 
When p > Pc the system is in a condensed phase. Here, 
the distribution function decays algebraically as p{m) ^ 
with an additional peak for condensates centered at 



,-6 



(15) 



The peak height scales as p(mcon) ^ N ^(In A^) ^/■^ and 
the width as Am ^ N^^"^. So the number of condensates 
A^con - A^p(»7icon)ATOcon scales as 



(16) 



The effect of the PE is evident. It breaks a macroscopic 
condensate into the number A'con of mesoscopic conden- 
sates of size rrican- The total mass of condensates is still 
macroscopic, i.e., Afcon^JT-con A^. 

The PE also plays an important role in the critical 
scaling ed, p — pc, where the mass distribution is given 
by p(to) ^ f{in). The maximum value of m, denoted as 
mc{N), will scales with A^. We will derive the scaling 
law for the maximum mass by using the extreme value 
statistics 12|, |l3|. Let F{mc) be the probability that the 
maximum mass is rric. It is given by 

F(m,) = a-^Nf{m,) a'' ^ f{m') 

\ m'=0 / 

where a — X]m=o fi^) ^ normalization con- 
stant. For large N and rric, it can be written as 
F{mc) ^ Nf{mc)exp[—NQ{mc)/a] where Q{mc) — 
dm' f{m'). Hence, the maximum mass mc{N) can be 
determined from the relation (5 (mc) ~ 1/A^. The integral 
with /(m) in Eq. ([T2|) yields that 
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FIG. 1: (Color online) Mass distribution function p(m) mea- 
sured in one-dimensional rings [filled symbols] and complete 
graphs [open symbols] with b — 4 and p = 4. The solid lines 
represent the curves for p{m) obtained from Eq. (|13p . 



where r(s, z) 
function |14| . 

oo) ~ z''"-'- 



= t" *di is the incomplete gamma 
It has the limiting behavior of F(s, z ^■ 



e ^ and r(s, z — > 0) 



r(.) 



Vs. 



m 

,1-fe 



(l+b)/2 -\r. 



Consequently, the function Q{m) behaves as Q{m) 

1 f^V' 

2 \2pN J 

g(m) ~ m^-" for (Am^) < A^. Solving Q{mc) 
self-consistently, we obtain that 



/{2pN) (;y^2) ^ ^ ^j^^ 

1/A^ 



{NlnNy/^ for b < 3, 

^1/(6-1) foj. 5>3^ 



(18) 



This is contrasted to the ZRP without PE (A = 0), where 
^ A^i/(^-i) for all values of 6 > 2 [l3. The results 
of Eqs. and justify the approximation ignoring 
0{{m/Mf) term in vim, m'). 

We have performed a numerical analysis to confirm 
our analytic results in one-dimensional rings and com- 
plete graphs with unbiased hopping. The single-site 
mass distribution function p{m) has been measured from 
Monte Carlo simulations. We have also evaluated p{m) 
numerically by enumerating the analytic expressions in 
Eqs. (HI]) and The normahzation factor Z{N,M) 

can be evaluated from the recursion relation Z{N, M) — 
E^^ofi^)Z{N-l,M-m). 

We present numerical data for p(m) at 6 = 4 and p = 4 
in Fig.[TJ The plots for one-dimensional rings overlap per- 
fectly with those for complete graphs. This suggests that 
it is likely that the particle species distribution is indeed 
random even in one-dimensional rings. There is a little 
deviation between the simulation data and the analytic 
results drawn with solid lines. The deviation comes from 
the approximation that 0{{m/N)'^) term is neglected in 
V. Since m is limited by the cutoff mpE N^^"^, the 
approximation is expected to be better as A^ increases. 
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FIG. 2: (a) mcon and A^con in a condensed phase at 6 = 4 
and p = 4. The fitting curves explained in the text are drawn 
with soHd fines, (b) The maximum mass mc{N) at the critical 
points p =!/(& — 2) at several values of b. Solid fines represent 
the analytic prediction of Eq. (|18|) . 

Numerical data shows that the deviation is already neg- 
ligible at TV = 1600, which justifies the approximation. 

The distribution function p{m) in Fig. [T] displays a 
broad peak representing condensates in the large m re- 
gion. The typical size mcon of condensates is given by 
the peak position, and the number iVcon of condensates 
is given by A'con = Nil. where fi is the spectral weight 
of the peak. The spectral weight can be estimated as 
^ ~ X]m>m p(to) where rrimin is the location of the 
local minimum of p{m). The numerical data thus ob- 
tained at 6 = 4 and p = 4 are presented in Fig. ^a). 
In order to compare with the analytic result of Eqs. ((15]) 
and ([TB|) . we fitted mcon{N) and Ncon{N) to the functions 
ai(7Vln(iV/a2))i/2 and asiN/ \ii{N/ai))^/^ , respectively. 
Numerical data are in good agreement with the fitting 
curves drawn with solid lines for large N. 

We have also examined the scaling behavior of the 
maximum mass mdN) along the critical line p = pc — 
l/(&— 2). The analytic theory based on the extreme value 
statistics predicts the scaling law in Eq. (fT8)) . We have 
measured mdN) at several values of b and present them 
in Fig. mb) . The solid lines represent the analytic result 
of Eq. p^ . Numerical data are in agreement with the 
analytic prediction for large values of N. 

In summary, we have shown that the pair exclusion can 
break a macroscopic condensate into multiple mesoscopic 
condensates. The pair exclusion is a weak interaction in 
that a particle interacts with only a single one. Such a 
weak interaction is negligible in a fluid phase. However, 
it becomes important during the course of condensation 
and makes a macroscopic condensate fragile. The meso- 
scopic condensation occurs in broader cases. Consider 
a pair repulsion that a pair occupation is rejected with 
probability A < 1. Such a system is also described by the 
same product state with Eq. PT|) . The cutoff mpE scales 
with the same exponent, hence the system displays the 
mesoscopic condensation for any nonzero values of A. 

The pair exclusion type interaction can play an impor- 



tant role in a broad range of phenomena. It is manifested 
as the self-loop constraint in network theory, whose ef- 
fects have been neglected in most studies. Percolation is a 
macroscopic condensation of sites into a single cluster. It 
will be interesting to study how the pair exclusion type 
interaction affects the nature of the percolation transi- 
tion. In fact, recent studies on the Achlioptas process [ll] 
reveal that a slight modification in a cluster growth rule 
changes a continuous percolation transition into a dis- 
continuous one [l5l - [l7j |. Another interesting application 
is an opinion dynamics modeled by a voter model fl8| . 
It may be interesting to study the effect of hostility be- 
tween individuals on consensus formation. Hopefully our 
work will initiate further studies in those directions. 
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